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The use of indirect turbulence measurements for real-time estimation of parameters in a 
linear longitudinal dynamics model in atmospheric turbulence was studied. It is shown that 
measuring the atmospheric turbulence makes it possible to treat the turbulence as a 
measured explanatory variable in the parameter estimation problem. Commercial off-the- 
shelf sensors were researched and evaluated, then compared to air data booms. Sources of 
colored noise in the explanatory variables resulting from typical turbulence measurement 
techniques were identified and studied. A major source of colored noise in the explanatory 
variables was identified as frequency dependent upwash and time delay. The resulting 
upwash and time delay corrections were analyzed and compared to previous time shift 
dynamic modeling research. Simulation data as well as flight test data in atmospheric 
turbulence were used to verify the time delay behavior. Recommendations are given for 
follow on flight research and instrumentation. 
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I. Introduction 

In operational flight conditions, an aircraft can encounter atmospheric turbulence and wind shears of varying 
magnitudes at any time, and many of these atmospheric turbulence encounters can occur in icing clouds. Real-time 
identification of dynamic models in atmospheric turbulence is a requirement for the Icing Contamination Envelope 
Protection ICEPro 1 system. ICEPro performs its function by making real time estimates of the stability and control 
characteristics of an ice contaminated aircraft in order to provide envelope protection control cues and messages to 
flight deck displays. By effective use of the cueing on these displays, a pilot is able to mitigate a potentially 
hazardous icing encounter by avoiding flight conditions and configurations that can result in an upset condition. 

The three most prevalent methods for identifying dynamic models from flight data are output error, equation 
error, and filter error. It has been shown that output error and equation error methods yield biased estimates of 
dynamic model parameters in the presence of atmospheric turbulence 2-3 . The output error method assumes 
measurement noise exists but not state noise, whereas equation error assumes that state noise exists but not 
measurement noise. Filter error methods represent the most general approach, where the dynamic system model has 
both state and measurement noise 4 . If filter error methods are used 5 " 9 , discriminating signal from noise is 
implemented through weighting matrices that represent assumed measurement and process noise covariances. In this 
case, adjustment of one or more tuning parameters must be done in simulation. For real-time applications, the 
required tuning and computational expense of filter error methods makes them difficult to apply in practice. 

Simulation studies accomplished for this work showed that if the atmospheric turbulence could be measured 
accurately, then real-time parameter estimation developed in previous work 10 ’ 11 could be applied directly to obtain 
accurate real-time parameter estimates in atmospheric turbulence, by treating the turbulence as a measured input in 
the aircraft dynamic model. This motivated a search for flight hardware that could accurately measure the 
turbulence. 

To explore the feasibility of the successful atmospheric turbulence simulation findings, commercial of the shelf 
(COTS) atmospheric turbulence sensors were investigated and compared against airdata booms. This approach 
makes it possible to treat the turbulence as a measured explanatory variable in the parameter estimation problem, 
and enables the use of equation error in the frequency domain, which does not involve tuning parameters, has very 
low computational requirements, and can be done in real time. 

Following the investigation on COTS atmospheric turbulence sensors, several potential sources of colored noise 
in the explanatory variables are investigated. These include: frequency response of airflow angle vanes, wing 
aerodynamic response in atmospheric turbulence, and the structural response of the air data boom. Numerical 
experiments using simulated and flight test data are included to support the investigation of each colored noise 
source. The simulated data come mainly from a linear DHC-6 Twin Otter simulation. Noise sequences similar to 
those observed in real flight data analysis and modeling were added to simulated aircraft responses. This approach 
makes it possible to identify practical problems and solutions regarding the parameter estimation, because the true 
parameter values are known from the simulation. Flight test results from DHC-6 Twin Otter aircraft are then 
compared and contrasted to numerical experiments for each of the identified error sources. Following this analysis, 
conclusions are made and recommendations are given for follow-on flight research and instrumentation. In general, 
all of the tools used to generate the results shown in the paper are available in a MATLAB software package called 
SIDPAC (System IDentification Programs for AirCraft) 10, 12 . 
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II. Atmospheric Turbulence Sensors 


A. Background 

In general, atmospheric turbulence is measured indirectly, because atmospheric turbulence is imbedded in the 
true airspeed measurement. The wind velocity vector V w relative to the surface of the Earth is found from the vector 
subtraction: 


n-=v,-v e 


( 1 ) 


where V t is the true airspeed and V e is the velocity of the aircraft with respect to the Earth, where the Earth is 
assumed to be an inertial reference. Because V w is computed as the difference between two large quantities, it is 

essential that V t and V e be measured as accurately as possible in order minimize errors in the computed V w . The 
most obvious difficulty in measuring atmospheric turbulence using a sensor mounted on the aircraft is removing the 
aircraft velocity with respect to the Earth. Additional difficulties arise because the aircraft distorts the airflow around 
itself. Once the wind velocity vector is calculated it can be treated as an explanatory variable in the parameter 
estimation problem. Any inaccuracies in the measurements or calibrations will affect the calculated wind velocity 
vector and in turn the parameter estimates. 

A number of atmospheric turbulence sensors have been developed over the years. Although successful 
measurements are possible using a simple 3 -axis accelerometer system to determine the response of the aircraft to 
gusts of wind 13 , this can be subject to large errors due to the mass of the aircraft damping the response to 
atmospheric turbulence. Notess 14 coupled accelerometers with fixed wind vanes and a pitot-static aircraft true air 
speed (TAS) sensor, whereas Telford and Warner 15 improved on this system with the introduction of the gyro- 
stabilized platform. Although fixed and rotating vanes 16, 17 have been used to measure airflow relative to the aircraft; 
it is now more common to use a differential pressure method 18-20 to avoid the vibration and frequency response 
problems of long airdata booms. Further refinement of the hemispherical pressure sensing head used to sense flow 
angles is described by Wyngaard 21 and Wood 22 . A more advanced instrument is the BAT Probe 23-25 which uses a 
nine hole pressure port formation in conjunction with differential GPS and fast-response accelerometers to give 
wind vector measurements at up to 50 Hz. 


B. Theory 

The wind velocity vector is found from the vector subtraction in Eq. (1). The following matrices are used to 
transform the true airspeed from the flight path axis V t to the earth axis (V ?yv y, E ,V^) 26-28 . The multiplication must 
be performed in the exact order /?, a , #, yr . 
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The matrix transformation from the flight path axis true airspeed to the earth axis (North, East, and Down) true 
airspeed is: 


3 

American Institute of Aeronautics and Astronautics 



{Yn + VwN ) 


X 

1 {V eE +V wl E ) 

> = L(i//)L(0)L(0)L(cc)L(j3)< 

0 

(YeD+V wD ) 


0 


where 


V wN = wind speed north 
V wE = wind speed east 
V wD = wind speed down 

From Eq. (4), the wind velocity components are calculated: 
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Embedded in the wind velocity vector is the atmospheric turbulent air velocity. The effect of atmospheric 
turbulence at a steady flight condition can be visualized by plotting a time history of the wind vector components 
and noting the variation in each component. 

This theoretical development assumes that the requisite variables are measured without error. Practical 
implementation of this theory necessitates ground and inflight calibrations to account for boom misalignment, 
up wash, side wash, angular rate error, boom bending, and response lag, among others. In typical flight test 
applications, the Euler attitude angles <9, y/ are measured using an inertial navigation system, while the airflow 

angles are measured using an airdata boom. Any inaccuracies in either measurement system will directly affect the 
computed wind velocity components. Special emphasis is placed on the word “computed” because the wind velocity 
components are not directly measured but computed using the vector subtraction in Eq. (5). 

C. COTS sensors 

A search was conducted to identify candidate COTS atmospheric turbulence sensors. Several sensors were 
identified using state of the art LIDAR 29 technology with an airflow angle accuracy of 1.0 degrees. In comparison, 
flight test air data booms with traditional air data technology are calibrated to an accuracy of 0.1 degrees. The cutoff 
accuracy for this work was set at 0.1 degrees. Unfortunately, the accuracy of LIDAR based sensors requires an order 
of magnitude improvement to be useful for the present work. Several additional sensors were identified using 
traditional air data technology 30, 31 . The AIMMS-20 probe by Avantech Research 31 was chosen for further 
investigation because it provides the most accurate wind determination of all the candidate sensors. The AIMMS-20 
probe incorporates state-of-the-art technology including differential GPS and miniaturized components, and makes 
use of Kalman Filter signal processing algorithms. Table 1 summarizes the AIMMS-20 performance. 

A quantitative comparison between the AIMMS-20 probe and typical flight test air data and inertial 
instrumentation will yield similar performance specifications for the measurements they have in common. The 
advantage of the AIMMS-20 probe is that it is packaged as a COTS system. The disadvantage of the AIMMS-20 
probe and all of the other COTS products identified is that atmospheric turbulence is not measured directly but 
computed. Hence, the COTS solutions are no different than applying the theory to flight test hardware, such as an air 
data boom and an inertial navigation system with differential GPS. Furthermore, any errors that affect an air data 
boom and an inertial navigation system will also affect the COTS atmospheric turbulence sensors. Therefore, it was 
found that current COTS atmospheric turbulence sensors do not offer an advantage over modern flight test 
hardware. 
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Table 1 AIMMS-20 probe performance 


Parameter 

1 Sigma Error 

Longitudinal acceleration, A x , (g) 

0.003 

Lateral acceleration, A Y , (g) 

0.003 

Normal acceleration, A z , (g) 

0.05 

Angle of Attack, a, (deg) 

0.10 

Angle of Sideslip, p, (deg) 

0.10 

Pitch rate, q, (dps) 

0.01 

Roll rate, p, (dps) 

0.07 

Yaw rate, r, (dps) 

0.01 

Pitch attitude, 0, (deg) 

0.20 

Roll attitude, (p, (deg) 

0.20 

Yaw attitude, \j/, (deg) 

0.20 

GPS velocity down, V D , (mps) 

0.10 

GPS velocity north, V N , (mps) 

0.10 

GPS velocity east, V E , (mps) 

0.10 

True airspeed, V T , (mps) 

0.20 

Wind velocity down, W D , (mps) 

0.30 

Wind velocity north, W N , (mps) 

0.30 

Wind velocity east, W E , (mps) 

0.30 


III. Problem Statement 

The following results are representative of the effect of atmospheric turbulence on aircraft dynamic model 
parameter estimates. Data were collected using the NASA Langley AirSTAR T-2 subscale generic transport aircraft 
during low and high atmospheric turbulence conditions. The T-2 aircraft is a 5.5 percent dynamically-scaled generic 
twin-engine jet transport aircraft. It is outfitted with flight test instrumentation and air data booms with flow vanes 
on each wingtip. Figure 1 shows a photograph of the aircraft in flight. 



Figure 7. T-2 subscale jet transport aircraft 
Credit: NASA Langley Research Center 
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Figure 2 shows parameter estimation results from T-2 flight test data using equation-error in the frequency 
domain. Each marker represents a different run, and only the last estimate of each real-time parameter estimation 
run is shown. A total of nine runs are represented; six in high turbulence and three in low turbulence. 
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Figure 2. Effect of atmospheric turbulence on T-2 modeling results 


Determining exactly what characteristic or effect from atmospheric turbulence causes the bias and scatter in the 
parameter estimation results is an ongoing research topic. It is known that colored noise in the explanatory variables 
or model outputs causes increased uncertainty and bias 2,3,10 in the parameter estimates, similar to what is observed in 
Figure 2 during repeated flight test maneuvers in high atmospheric turbulence. It is unknown whether the turbulence 
causes colored noise in the explanatory variables and/or outputs, or if the observed bias and uncertainty in the 
parameter estimates in turbulence are from some other cause. 

A. Hypothesis 

Colored noise in the explanatory variables is the source of bias and increased uncertainty in the parameters 
estimated in high atmospheric turbulence. 

B. Sources of colored noise 

Possible sources of colored noise in the explanatory variables include: frequency response of airflow angle 
vanes, frequency response of the boom, wing aerodynamic response to atmospheric turbulence, pitot-static 
frequency response, and airframe structural dynamics, among others. Since the airflow angles a and p are the 
explanatory variables with the largest contribution in the dynamic modeling process, the focus will be on these 
quantities. The effects of the frequency response of airflow angles vanes, the wing aerodynamic response in 
atmospheric turbulence, and the structural response of the air data boom will be considered using simulation and 
flight test data. 
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1. Frequency response of airflow angle vanes 

In order to investigate the frequency response of airflow angle vanes, an airflow vane sensor model was studied. 
A natural frequency of co n = 9 Hz and a damping ratio <^ = 0.35 were chosen as representative values for a 
dynamic model of the airflow vanes 32 ' 34 . Typical natural frequency and damping ratio values range from 5-20 Hz 
and 0.2-0. 6 respectively. Eq. (6) represents a second order continous-time transfer function. 


G(s) 


CO„ 


s 2 + 2 dfco n s + CO 2 


( 6 ) 


The transfer function representation does not support non-zero initial conditions; therefore, the vane sensor model 
was converted to a state space representation which supported a trim a initial condition. 

2. Wing aerodynamic response in atmospheric turbulence 

A DHC-6 Series 300M Twin Otter research aircraft was used to study the aerodynamic effects on atmospheric 
turbulence measurements 35 . Details of the aircraft and its instrumentation system are found in Ref. [35]. A flush air 
data system mounted on the nose of the Twin Otter was used as the primary source of airdata measurements. Flush 
air data systems avoid the vibration and frequency response problems of long booms mounted on the nose of the 
research aircraft but require extensive calibrations 36 . 

It was determined 35 that a major flow distortion, which seriously affects the vertical wind velocity measurements 
near the nose of an aircraft, is the frequency dependent upwash and time delay. These are shown to be the result of 
not taking into account the decrease in upwash correction with eddy frequency caused by the response of the wing 
vortex system to atmospheric turbulence, and the longitudinal distance between the angle of attack sensor and the 
leading edge of the wing, respectively. 

The unsteady aerodynamic analysis detailed in Ref. [35] with flight data from NASA’s DHC-6 Series 100 Twin 
Otter aircraft was used to calculate the angle of attack upwash and time delay as a function of frequency. The 
required upwash correction and time delay is summarized in Figure 3 and Figure 4 during flight in light and 
moderate turbulence. Similar results are found in Ref. [35]. 
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Figure 3 . Angle of attack upwash correction versus frequency 

Figure 3 represents the estimated angle of attack scale factor correction vs. frequency. The frequency scale is the 
reduced frequency k b based on the half wing span where co is the angular frequency, b is the wingspan, and U is the 

airspeed. The angle of attack upwash correction (l + k u ) 1 is explained by examining a sensor model: 
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(7) 


a M = {\ + k u )a T +a B 


where 


a B = bias angle of attack 
a M = measured angle of attack 
a T = true angle of attack 


In Figure 3, the vertical axis (l + k u ) Represents the correction required to calculate the true angle of attack from 
measured data. 


a T=( l + k u) \ a M- a B) ( 8 ) 

At low frequencies (representing flight in a calm atmosphere) the scale factor is approximately 0.75 and 
increases to 0.90 as k b approaches 2.0. A k b value of 2.0 (using the DHC-6 Twin Otter wing span and 110 knot 
cruise airspeed) corresponds to approximately 2.0 Hz. Therefore, over the frequency range of interest for rigid body 
aircraft dynamics (0.1 to 2.0 Hz), the scale factor changes from 0.75 to 0.90 due to atmospheric turbulence. 
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Figure 4. Time delay versus frequency 

Figure 4 depicts the time delay as a function of frequency. The time delay is caused by the separation between 
the nose boom mounted airflow vanes (measurement position) and the wing tip (lift and up wash generation), along 
which turbulence eddies travel. The magnitidue of the time delay represents the correction required to calculate the 
true angle of attack from measured data. In the frequency range of interest for rigid body aircraft dynamics (0.1 to 
2.0 Hz), the time delay is approximately 0.1 seconds. A general description of unsteady aerodynamic theory at the 
nose of the aircraft and its calculations are found in the seminal paper 35 . 

A data compatibility 10 analysis was used to reconstruct the angle of attack from flight test measurements and 
compare it to the measured angle of attack. Three runs were evaluated at similar airspeeds (120 knots) and aircraft 
configurations. The first run represents a calm atmosphere; the next two runs represent atmospheric turbulence. 
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Time Delay in Angle of Attack Regressor Run 1 Calm Atmosphere 



Figure 5. Angle of attack time delay in atmospheric turbulence 


In Figure 5, run 1 in calm atmospheric conditions shows a time delay of approximately 0.03 seconds. Runs 2 and 
3 in atmospheric turbulence show a time delay of approximately 0.13 seconds. This agrees with the fact that the time 
it takes to travel 21 feet (longitudinal distance from the airflow angle measurement position to the center of gravity) 
at 120 knots is approximately 0.1 seconds. The difference is attributed to forced (elevator doublets) oscillations in 
run 1 and a combination of forced and atmospheric turbulence oscillations in runs 2 and 3. This increase in time 
delay in atmospheric turbulence is a source of degradation in the stability and control parameter estimates. 

Another problem with using inertial measurements to reconstruct data such as angle of attack is that the 
reconstruction assumes the atmosphere is fixed relative to the earth axes. When there are gusts or winds, this 
assumption is not valid, and the reconstructed data have errors. Since the airflow angles are measured at the nose of 
the aircraft and the inertial measurements are taken at the center of gravity it is expected that the measured angle of 
attack will lead the reconstructed inertial angle of attack. The reconstructed angle of attack leads the measured angle 
of attack and is in error in Figure 5. This is attributed to not accounting for gusts and winds in the data compatibility 
analysis. 

3. Structural response of air data boom 

In order to investigate the inflight structural response of the airdata boom, airflow angle flight data collected in a 
DHC-6 Twin Otter during multiple parameter etstimation maneuvers was analyzed in the frequecy domain using 
Fourier sine series coefficients. 

Figure 6 shows the frequency response data for the entire frequency range, up to the Nyquist frequency at 25 Hz. 
The longitudinal and lateral natural frequencies for the dominant structural mode of the air data boom are 
approximately 8 and 9 Hz. These natural frequencies are associated with the angle of attack and angle of sideslip 
measurements respectively. In Figure 6, the Fourier coeffcients near 8 and 9 Hz are at the noise floor. Therefore, the 
structural response of the nose boom on which the angle of attack sensor is mounted and from fuselage bending is 
not a source of colored noise in the angle of attack. Similar results were calculated for angle of sideslip but the 
results are not shown. 
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Angle of Attack 



Figure 6. Fourier sine series coefficients for measured angle of attack 


IV. Real-Time Parameter Estimation 

Time-domain methods are generally not good for real-time applications due to the iterative nature of the 
solutions, or trial and error weighting of the formulations, or the fact that estimated parameter error bounds are 
either inaccurate or require rather involved correction calculations 10 . However, if equation error in the frequency 
domain is used, the formulation accounts for state noise and is also robust to measurement noise, in addition to 
providing accurate error bounds without correction calculations. This is due to the automatic filtering inherent in 
using a limited bandwidth for the recursive Fourier transform 11 since the rigid-body dynamics of typical piloted 
aircraft lie in a relatively narrow frequency band of approximately 0.01-2.0 Hz. 

The real-time parameter estimation algorithm being employed was developed by Morelli 11 , and is based on a 
frequency domain approach implemented in version 2.0 of System IDentification Programs for AirCraft 
(SIDPAC). 10, 12 The approach uses a linearized model of the vehicle dynamics, where A, B , C, and D are matrices 
containing stability and control derivatives. 


x(t} - Ax(t^ + Bu(t} x(0) = 0 

y (V) = Cx(t) + Du(t) 


( 9 ) 


In Eqs. (9), the stability and control derivatives to be estimated can be dimensional or non-dimensional 10 . Time- 
varying estimates of the stability and control derivatives are determined using measured aircraft control surface 
deflections and measured states from air flow angle data and angular rate data. Outputs are aircraft states and 
translational accelerations. To determine the model parameter estimates, a cost function is formulated in the 
frequency domain from the Fourier transform of the model in Eq. (9), 


joj.x(oj) = Ax(oj) + Bit (fo>) 
y(co) = Cx(rv) + Du 


(10) 
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The cost function is 
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for m frequencies 


( 11 ) 


where ^ and b k are the kth row of matrices A and respectively, and x k ( n ) = x k ( co n ), (zz) = zz^ (co n ) . Each 

line in Eq. (10) can be analyzed separately in this way, which implements an equation-error formulation 11 . The least 
squares cost function in each case can be formulated as 


J A(Y-X9) f (Y-X6) 


where 
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( 12 ) 


(13) 


and the unknown parameters from A and B are contained in the parameter vector 9. The least squares parameter 
vector estimate is obtained as the value of #that minimizes the cost function in Eq. (12). 10 


9 = 


1 Re[x f Y^ 

cov[9) = E (<9-6»)(<9-6») T 


Re 




(14) 


Performing the identification in the frequency domain is advantageous for this application because of 
computational efficiency, robustness to noise and data dropouts, and reliable confidence bounds from the covariance 
matrix. The required Fourier transforms are computed using a recursive formulation, and the calculations in Eq. (14) 
are done periodically (typically at 1 or 2 Hz) to implement real-time parameter estimation, see Refs. [10] and [11] 
for details. 


V. Aircraft 

The DHC-6 Series 100 Twin Otter is a Canadian 19 -passenger Short Takeoff and Landing utility aircraft 
developed by de Havilland Canada. Figure 7 shows a photograph of the aircraft in flight. 

The aircraft has two 550 SHP Pratt and Whitney PT6A-20A turbine engines driving three bladed Hartzell 
constant speed propellers. The flight controls are mechanically operated through a system of cables and pulleys. 
Control surfaces consist of elevator, ailerons, rudder, and wing flaps. Table 2 contains mass and geometry properties 
of the Twin Otter aircraft. 
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Credit: NASA Glenn Research Center 


Table 2. Physical properties of the Twin Otter aircraft 


Property 

Value 

m, slug 

280 

Inertia 


iy, SlUg-fd 

16,000 

Wing 


Area, ft 2 

422.5 

Span, ft 

65 

Mean Aerodynamic Chord, 
ft 

6.5 


A. Aircraft Simulation 

The flat earth, body axes, longitudinal three degree of freedom model equations can be classified as the force 
(15), moment (16), kinematic (17), and navigation (18) equations. This formulation was chosen to enable both batch 
simulations and pilot in the loop simulations with realistic out-the- window visuals. The dynamics were reduced to 


the three longitudinal degrees of freedom by setting lateral and directional values 
order to simplify the analysis. 

v = p = r = (j>=r = p = O 10 in 

u = qw— g sin 6 + (A + T x ) / m 

(15a) 

w — qu-g cos 0 + (Z + T z )l m 

(15b) 

fq — M + 

(16) 
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( 17 ) 


h = usinO-wcosd (18) 

Nondimensional aerodynamic force and moment coefficients for an aircraft can be computed from flight 
measurements as follows 10 : 

C x — (ma x —T x )/qS (19a) 

C z =(ma z -T z )/qS (19b) 

C M =[I Y q-M T ]/qSc (19c) 

For local modeling over a short time period, the force and moment coefficients computed from Eqs. (19) can be 
modeled using linear expansions in the aircraft states and controls: 


Cx C Xa Aa + C Xq +C XSe A8e + C Xc 


Cz ~ C Za Act + C Zq ^ + C ZSe A5e + C Zo 


C M =C Ma Aa + C Mq -^- + C M5e ASe + C Mo (20c) 

The A notation indicates perturbation from a reference condition. In Eq. (20c), C mo represents the 

nondimensional pitching moment at a reference condition, and similarly for the other expansions. For short periods 
of time, the stability and control derivatives are considered as constant model parameters to be estimated from flight 
data. 

B. Atmospheric Turbulence model 

Atmospheric turbulence was tested as calm, light, moderate, and severe. Atmospheric turbulence levels and 
implementation were defined per MIL-F-8785C and were based on the Dryden turbulence spectrum. The 
atmospheric turbulence velocity components were generated randomly within the guidelines of the Dryden spectrum 
using individual realizations of random number sequences and added to the aircraft velocity components to generate 
the total velocity components. This resulted in unique sequences of random numbers for each test. The three 
velocity components from the atmospheric turbulence calculations were added to their appropriate body frame 
velocity components (w, v, and w) calculated from the aircraft dynamic equations. 

u ~ u aerodynamics u wind u turbulence (2 1 n) 

v ~ V aerodynamics v wind T turbulence (21b) 

^ — W aerodynamics ^ wind ^ turbulence (2 1 c) 
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Medium/High Altitude Turbulence Intensities (Probability of Exceedance) 



Figure £. 37 MIL-F-8785c atmospheric turbulence velocities 

The root-mean- square turbulent wind speeds are also shown Figure 8 as a function of altitude. According to 
MIL-F-8785C, the probability of exceeding light turbulence levels is between 10" 1 and 10" 2 , for moderate turbulence, 
the probability is approximately 10‘ 3 , and for severe turbulence it is approximately 10‘ 5 . Therefore, severe turbulence 
levels are not likely to be encountered. 


VI. Simulation Results 

The hypothesis was tested using a simulation of the DHC-6 Series 100 Twin Otter and the following test conditions: 

1) Baseline runs with perfect angle of attack measurements 

2) Runs with a representative frequency response model for the airflow angle vanes 

3) Repeat test condition 2) with varying angle of attack magnitude and time delay as a function of frequency 
The following parameters are common to all three test conditions: 

- Aircraft trimmed at 2,500 feet and 110 knots 

- 20 second runs 

- Real-time parameter estimation method 10, 11 , described in Section IV 

- Orthogonal multi-sine inputs 10 invoked from 1 to 11 seconds 

- Selected final model parameter estimates from each 20-second run 

- 200 runs per atmospheric turbulence intensity (calm, light, moderate, and severe) 

- Fourier transform frequency spacing 0.04 Hz 

- Fourier transform frequency range [0.1 2.0] Hz 

- Signal to noise ratio of 50 

- 50 Hz data rate 
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The following comments are common to all of the figures associated with test conditions one through three: 

- Green line- truth value of the parameter 

- Red dashed line - confidence interval using ±10% of truth value of the parameter 

- Red diamond - parameter estimate 

- Blue bars - estimated two sigma error bounds 

In the dynamic modeling process, no attempt was made at estimating the axial force coefficients, because the aircraft 
excitations did not appreciably change airspeed, resulting in a low signal-to-noise ratio for the axial force coefficient 
parameter estimation. In total, six coefficients were estimated, C Za ,C Zq ,C ZSe ,C Ma ,C Mq , and C MSe . 

A. Test Condition 1 

As show in Figure 9, the baseline model estimates for C Za resulted in improved estimates with increasing 

atmospheric turbulence. This is not representative of actual flight data and is attributed to perfect a measurement 
given the additional excitation caused by atmospheric turbulence. 
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Figure 9. Test condition 1 - effect on C Za 
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As shown in Figure 10, the baseline model estimates for C MSe resulted in greater scatter as turbulence intensity 
increased. This is caused by noisy explanatory variables 38 from the atmospheric turbulence implementation in 
Eq. (21). Similar patterns were observed with the remaining coefficients C Zq ,C ZSe ,C Ma , and C Mq , although the 

results are not shown. 
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Figure 10. Test condition 1 - effect on C MSe 
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B. Test Condition 2 

Estimates for C Za based on the vane sensor model ( co n - 9 Hz and — 0.35 ) caused degradations and biasing 
in the estimates with increasing turbulence intensities, as shown in Figure 11. This is representative of actual flight 
data. 
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Figure 11. Test condition 2 - effect on C Za 
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As shown in Figure 12, the baseline model estimates for C MSe resulted in greater scatter as atmospheric 

turbulence intensity increased. This is caused by noisy explanatory variables 38 from the atmospheric turbulence 
implementation in Eq. (20), and the vane sensor model. Similar patterns were observed with the remaining 
coefficients C Zq ,C ZSe ,C Ma , and C Mq . Results are not shown. 


Calm Atmosphere Light T urbulence 




Figure 12. Test condition 2 - effect on C MSe 

C. Test Condition 2 Expanded Testing 

Although no quantitative analysis was done, qualitatively, small differences were observed between a perfect 
and typical vane frequency response in test conditions 1 and 2. The vane dynamics do not substantially contribute to 
the scatter in the estimates, because all of the added dynamics are outside the range of frequencies analyzed. 

To further validate this claim, additional testing was carried out using the same simulation. A total of 121 
different combinations of flow angle vane natural frequency and damping ratios were tested. The natural frequency 
values tested were 5 through 15 in 1 Hz increments. The damping ratio values tested were 0.01, 0.05, 0.1, 0.15, 0.2, 
0.25, 0.3, 0.35, 0.4, 0.45, and 0.6. Since the AIMMS-20 probe approaches a gain of one and zero phase shift, it is 
characterized by the results of test condition 1. The AIMMS-20 probe offers superior frequency response 
performance to mechanical vanes by using 1) a 5 hole pressure probe , 2) short tube lengths (air data tubes are 
typically less than 5 inches long), and 3) the resonance frequency of the silicon pressure diaphragm is in the 
kilohertz region. The following figures depict the results of 121 combinations, times 4 turbulence levels, times 200 
runs per case. This totaled 96,800 runs. Figure 13 depicts the results of 96,800 runs forC Za . As an example, the 
bottom left corner block in severe turbulence is light blue. The color represents the number of exceedances of 
C Za from ±10% of truth. The remaining 120 blocks are dark blue representing less than 10 exceedences. Figure 13 
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Vane Natural Frequency Vane Natural Frequency 


suggests that any combination of natural frequency and damping ratio tested provides adequate performance. Similar 
results were obtained for C Ma and C Mq . Results are not shown. 
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Figure 13. Number of exceedances for C Za 
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Figure 14 shows that independent of the natural frequency and damping ratio tested, accurate estimates of C z§e 
are not feasible in atmospheric turbulence. This is expected since the normal force produced by elevator deflection , 
C z5e , is an order of magnitude less than the normal force produced by angle of attack, C Za . 
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Figure 14. Number of exceedances for C Mq 
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Figures 15 shows that independent of the natural frequency and damping ratio tested, accurate estimates of 
C MSe are feasible in atmospheric turbulence. This is expected due to the elevator control authority. 
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Results for C Zq , not shown, demonstrate that independent of the natural frequency and damping ratio tested, 

estimates are not feasible in any level of atmospheric turbulence. For the remaining stability and control estimates 
the results show that most combinations of natural frequency and damping ratio tested provide adequate 
performance in calm, light, and moderate turbulence, and severe atmospheric turbulence. 
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D. Test Condition 3 

As shown in Figure 16, estimates for C Za resulted in degraded estimates with increasing atmospheric turbulence 
for varying angle of attack magnitude and time delay as a function of frequency. The degradation is seen as biasing 
in the estimates and increased scatter and uncertainty as the error bars increase. This is representative of actual flight 
data. 
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Figure 16. Test condition 3 - effect on C Za 
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As shown in Figure 17, the estimates for C z§e also resulted in degraded estimates with increasing atmospheric 
turbulence. The degradation is seen as biasing in the estimates and increased scatter and uncertainty as the error bars 
increase. The increase in biasing and uncertainty compared to C Za estimates is caused by the relatively small 
contribution of the elevator term to the normal force coefficient. 
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Figure 17. Test condition 3 - effect on C z5e 
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Figure 18 shows that the estimates for C Ma degraded with increasing atmospheric turbulence. The degradation 

is seen as biasing in the estimates and increased scatter and uncertainty as the error bars increase. This is 
representative of actual flight data. 
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Figure 18. Test condition 3 - effect on C Ma 
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As shown in Figure 19, the estimates for C Mq degraded with increasing atmospheric turbulence. The 

degradation is seen as biasing in the estimates and increased scatter and uncertainty as the error bars increase. The 
estimates for C Mq resulted in the largest biasing of all the aerodynamic estimates considered because angular 

velocity model parameters are most susceptible to time delays 39 . 
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Figure 19. Test condition 3 - effect on C Mq 
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Figure 20 shows that the estimates for C MSe degraded with increasing atmospheric turbulence. The degradation 

is seen as a scatter in the estimates and increased scatter and uncertainty as the error bars increase. This is 
representative of actual flight data and correlates with the flight results presented in Figure 2. 




Figure 20. Test condition 3 - effect on C MSe 

Corrupting the angle of attack measurement by a frequency dependent magnitude and time delay affected all of 
the estimated stability and control coefficients. The estimates of C Ma and C MSe exhibited moderate biasing and 
increased scatter and uncertainty, C Za and C Mq estimates had large biases and increased scatter and uncertainty, 

and C z§e estimates showed the largest biases and increases in scatter and uncertainty. Previous work on the effects 
of time-shifted data on flight determined stability and control estimates 39 using the output-error method in the time 
domain demonstrated identical trends in the biasing and increased uncertainty due to atmospheric turbulence. 
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VII. Conclusions 

Several practical issues related to atmospheric turbulence in real-time estimation of parameters in a linear 
longitudinal dynamics model were examined and discussed. These issues included the use of COTS sensors and 
identifying sources of colored noise in the explanatory variables. It was determined that current COTS sensors 
compute the atmospheric turbulence and that any errors affecting an airdata boom and inertial navigation system 
will most likely affect COTS sensors. Therefore, current COTS atmospheric turbulence sensors do not offer an 
advantage over typical flight test hardware. The frequency response of the airflow vanes, wing response to 
atmospheric turbulence, and the structural response of the air data boom were identified as sources of colored noise. 
The frequency response of the airflow vanes for all combinations of natural frequency and damping ratio were 
adequate except for C Zq and C z§e estimates in atmospheric turbulence. This is attributed to noisy explanatory 

variables in the atmospheric turbulence implementation and relatively low aerodynamic contributions as compared 
to C Za . Therefore, the airflow vane dynamics themselves do not substantially contribute to the scatter in the 
estimates, because all of the added dynamics are outside the range of frequencies analyzed. A major source of 
colored noise is the frequency dependent upwash and time delay induced by the wing-bound vortex system and the 
longitudinal separation between the angle of attack measurement and the wing. The time delay in the angle of attack 
is the greatest contributor to biasing and increased scatter and uncertainty in the parameter estimates. 

This approach enables the use of equation error methods in the frequency domain, without the need for tuning 
parameters and complex calculations typical of filter error methods. This is accomplished by accounting for the 
sources of the colored noise in the measurements instead of treating the source(s) as unknowns and accounting for 
them through the use of process noise. 

Practical issues were examined using data from a Twin Otter DHC-6 longitudinal linear simulation, with realistic 
noise sequences added to the computed aircraft responses. This allowed a clear view of the effect of each source of 
colored noise to the modeling problem, because the true values of the model parameters were known. Flight test data 
from the NASA Glenn DHC-6 Twin Otter aircraft was used to validate the effect of the identified colored noise 
sources. 

Based on these findings, several practical recommendations for atmospheric turbulence flight test are suggested 
as well as areas for future study. First, mount the airdata boom at the nose of the aircraft. If the airdata boom must be 
mounted on the wingtip, additional analytical techniques are needed to separate the aero elastic effect of the wing 
from the low frequency aircraft response. Second, if the airdata boom is mounted on the nose of the aircraft, an 
upwash and time delay calibration as a function of frequency is required. Third, for frequencies of interest, 0.1 to 2.0 
Hz, either airflow vanes or a pressured based measurements are adequate. Suggested areas for future study include 
the airspeed frequency response, model structure determination, lateral/directional real-time parameter estimation, 
as well as applying the frequency dependant upwash and time delay to obtain calm air estimates in atmospheric 
turbulence. 
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